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摘要 : 由 于 特 高 压 直 流 输电 线 周 围 存在 空间 电荷 ， 直 流 输电 线 附 近 电 场 分 布 复杂 ， 
对 附近 的 绝缘 子 办 络 存 在 影响 。 本 文 基于 上 流 有 限 元 方法 计算 模拟 了 500kV 单 极 高 压 
直流 线路 附近 的 离子 流 场 ， 采 用 友 代 收敛 控制 算法 ， 以 保证 迭代 收 你 ， 研究 了 空间 电 
荷 的 存在 对 500kV 直流 绝缘 子 沿 面 参数 的 影响 。 
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Abstract: The space charge around ultra high voltage dc transmission line makes 
the electric field distribution nearby complicated, and has an effect on the insulator 
flashover nearby. The ion flow field near the S00kV single-electrode high voltage dc line 


is simulated based on the upper-class finite element method in this paper, and an iterative 


convergence control algorithm is adopted to ensure iterative convergence. The influence of 


space charge on the surface electric field parameters of 500 kV dc insulators is studied. 
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高 压 直 流 输电 技术 在 远 距 离 、 大 容量 输电 方面 
有 较 大 的 优势 "在 我 国电 网 中 发 挥 着 重要 的 作 
用 。 然 而 ， 高 压 直 流 输电 的 发 展 却 面 临 较 大 的 技术 
难题 ， 与 交流 输电 线路 相 比 ， 高 压 直 流 输电 线路 周 
围 存 在 的 空间 电荷 导致 线路 附近 电场 强度 与 交流 线 
路 有 明显 区 别 ， 对 附近 的 电器 设备 的 电场 分 布 也 有 
较 大 的 影响 ， 由 于 绝缘 子 距离 输电 线 十 分 接近 ， 受 
影响 较 明 显 ， 且 绝缘 子 表面 电场 电位 数值 对 电网 稳 
定 运行 具有 十 分 重要 的 现实 意义 。 因 此 ， 有 必要 对 
直流 线路 附近 的 离子 流 场 进行 研究 。 

同时 ， 直 流 输电 线路 运行 时 ， 如 果 导 线 表 面 的 
电场 强度 超过 大 气 的 击 穿 场 强 ， 导 线 表 面 会 发 生 电 
量 放 电钻 ， 产生 空间 电荷 ， 线 路 附近 的 电荷 会 在 空 
间 产 生 电 场 ， 使 直流 线路 和 交流 线路 具有 不 同 的 电 
场 分 布 。 高 压 直 流 线 路 在 运行 时 ， 会 发 生 电 坚 放电 
现象 ， 而 且 直流 输电 线路 与 交流 线路 在 发 生 电 晕 现 
象 时 ， 附 近 的 空间 电荷 分 布 存在 很 大 区 别 。 由 于 交 
流 输电 线 的 线路 电压 随时 间作 正弦 变化 ， 因 此 当 交 
流 线 路 产生 电 吧 时， 在 上 半 周 期 产生 的 空间 电荷， 
在 下 半 周 期 时 ， 由 于 线路 极 性 发 生 改变 ， 又 被 导线 
吸引 ， 因 此 ， 交 流 输 电线 附近 的 空间 电荷 仅 在 导线 
附近 做 往返 运动 ， 在 各 相 输 电线 线 之 间 和 相 导 线 与 
大 地 之 间 的 范围 内 基本 不 存在 空间 电荷 。 而 直流 输 
电线 路 发 生 电 晕 时 ， 由 于 输电 线 的 电压 极 性 恒定 ， 
在 极 导线 发 生 电 尝 产生 的 带电 离子 中 ， 与 导线 极 性 
相反 的 离子 向 导线 移动 ， 而 与 导线 极 性 相同 的 离子 
将 远离 导线 ， 沿 所 受 电场 力 方向 移动 ， 因 此 ， 在 各 
相 输 电线 线 之 间 和 相 导 线 与 大 地 之 间 的 范围 内 会 存 
在 带电 离子 。 
在 整个 输电 线路 上 ， 绝 缘 子 是 其 中 应 用 广泛 的 
装备 之 一 ， 电 网 的 安全 与 稳定 运行 与 绝缘 子 运行 状 
况 息息相关 。 由 于 直流 输电 线路 存在 离子 流 场 问 题 ， 
绝缘 子 串 的 周围 存在 着 空间 电荷 ， 使 其 电压 及 电场 
强度 分 布 特性 与 交流 绝缘 子 发 生变 化 。 为 了 保障 高 
压 直 流 线 路 绝缘 子 在 线路 正常 运行 ， 减 少 内 络 事故 ， 
在 设计 时 ， 需 要 分 析 绝 缘 子 的 沿 面 参数 。 由 于 直流 
线路 存在 离子 流 场 问题 ， 这 将 导致 绝缘 子 串 的 周转 
(尤其 是 高 压 端 ) 存在 着 大 量 的 空间 电荷 ， 这 些 空 间 
电荷 使 绝缘 子 的 沿 面 电 压 及 电场 强度 分 布 特性 与 交 
流 绝缘 子 不 同 。 当 绝缘 子 沿 面 电 压 分 布 不 均匀 的 情 
况 严重 时 ， 绝 缘 子 处 于 高 场 强 一 端的 部 件 容易 发 生 
老化 ， 绝 缘 能 力 受 到 破坏 ， 当 局 部 发 生 电 尝 时 ， 绝 
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缘 子 的 绝缘 能 力 受到 影响 更 大 ， 这 将 严重 影响 电网 
安全 运行 。 为 使 直流 输电 线路 更 安全 地 运行 ， 需 要 
对 直流 绝缘 子 沿 串 电压 、 电 场 分 布 计算 方法 进行 研 
究 。 这 种 研究 可 以 更 深入 地 了 解 直流 线路 中 绝缘 子 
的 性 能 ， 有 助 于 更 好 地 针对 直流 绝缘 子 的 特点 进行 
优化 设计 ， 同 时 ， 对 于 进一步 研究 绝缘 子 的 沿 面 放 
电 和 污 闪 提供 了 支持 ， 对 于 直流 线路 绝缘 子 的 优化 
设计 和 稳定 运行 具有 重要 意义 ， 通 过 这 些 优化 ， 可 
以 逐步 确保 直流 输电 线 的 安全 运行 ， 从 而 最 终 解 决 
高 压 直流 输电 的 技术 难题 。 

本 文 对 描述 直流 输电 线路 离子 流 场 的 Poisson 
方程 和 电流 连续 性 方程 进行 了 求解 ， 提 出 了 改进 的 
迭代 方法 ， 基 于 Kapzov 假设 ,研究 了 直流 输电 线 
的 离子 流 场 及 风速 对 其 影响 ， 同 时 使 用 Matlab 软 
件 进 行 编程 ， 采 用 有 限 元 法 ， 计 算 了 空间 电荷 对 
500kV 绝缘 子 表面 电场 和 电位 分 布 的 影响 。 


2 ”数学 模型 


2.1 离子 流 场 控 制 方程 

特 高 直流 输电 线路 在 工作 时 ， 由 于 发 生 电 尝 现 
象 ， 附 近 存 在 大 量 的 带电 离子 。 带 电离 子 的 运动 形 
成 离子 流 。 考 虑 风速 的 影响 ， 由 高 压 直 流 输电 线 的 
具体 运行 环境 及 麦克 斯 韦 方程 组 可 以 得 出 离子 流 场 
的 控制 方程 如 下 : 


Vig=-(p, -Pp_)/eo (1) 

E=-Vg (2) 
7,=p:(K.E+w) (3) 
7-=p-(K_E-w) (4) 
VJ,=-Rp.p-/e (5) 
VJ =-Rp,p_/e (6) 


式 中 , g 为 电位 ; p;、p- 分 别 为 正 、 负 空间 电荷 窗 
度 ; so 为 空气 介 电 常数 ， 为 合成 电场 强度 ，J,、J 
分 别 为 正 、 负 离子 电流 密度 ; K,、K_ 分 别 为 正 、 负 
离子 迁移 率 ; w 为 风速 ，R 为 离子 复合 系数 ; e 为 电 
子 电量 。 
2.2 基本 假设 

为 更 加 方便 地 对 离子 流 场 进行 分 析 ， 本 章 在 分 
析 中 采用 如 下 假设 : 

(1) 高 压 直 流 输电 线 周 围 空间 充满 了 两 种 极 性 
的 电荷 ， 由 于 电离 层 厚 度 远 远 小 于 导线 对 地 高 度 和 
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极 间距 离 ， 因 此 对 其 厚度 忽略 不 计 。 

(2) 计算 时 认为 导线 电 尝 已 经 稳定 ， 不 考虑 其 
中 的 暂 态 变化 过 程 。 

(3) 计算 时 导线 附近 的 电 尝 保持 稳定 ， 起 举 场 
强 不 发 生变 化 ， 即 Kaptzov 假设 。 

(4) 离子 迁移 率 上 和 复合 系数 了 与 其 他 变量 
关 ， 作 为 常数 。 

(5) 不 考虑 空间 电荷 扩散 ， 尽 管 空间 电荷 存在 
扩散 现象 ， 但 是 与 其 所 受 电场 力 相 比 ， 基 本 可 以 名 
略 。 因 此 忽略 电荷 扩散 不 产生 过 大 误差 。 

(6) 在 加 入 风速 条 件 时 ， 忽 略 风速 变化 过 程 ， 
风速 为 常量 且 稳 定 不 变 。 

(7) 不 考虑 杆 塔 影响 和 电 尝 分 布 的 不 均匀 性 ， 
使 三 维 问题 变 为 二 维 问题 。 

在 本 文 的 计算 中 ， 电 晤 起 始 场 的 大 小 通过 式 
(1) 计算 获得 地 面 合成 电场 的 约束 方程 为 


到-29sodl 


a 


式 中 , 了 为 导线 半径 ; 疡 为 导线 表面 粗糙 系数 ，6 为 
空气 的 相对 密度 。 


3 ”模型 求解 


在 求解 直流 输电 线路 离子 流 场 时 ， 需 要 利用 
用 边界 条 件 以 求解 泊 松 方程 和 电流 连续 性 方程 式 
(1) ~ 式 (6)。 其 中 ， 由 于 每 个 节点 的 空间 电荷 密 
度 p:、p- 和 电位 p 均 为 待 求 量 ， 本 文采 用 迭代 方法 
进行 求解 。 在 计算 时 ， 首 先 根据 同 轴 圆柱 的 电荷 密 
度 求 解 公式 给 出 各 节点 的 近似 电荷 密度 初 值 ， 根 据 
所 给 的 电 痊 密度 和 边界 条 件 计算 出 各 节点 的 电位 及 
场 强 ， 之 后 再 由 各 点 的 电场 强度 ， 求 解 电流 连续 性 
方程 ， 得 到 所 有 市 点 的 空间 电 奏 密度 ， 再 根据 新 的 
空间 电 痊 密度 分 布 求解 新 的 电位 及 其 场 强 ， 直 到 电 
场 强度 和 所 有 节点 的 空间 电荷 密度 的 计算 结果 同时 
满足 给 出 的 收敛 条 件 ， 则 这 个 计算 结果 就 为 所 求 数 
值 解 。 程 序 计算 流程 如 图 1 所 示 。 

本 文 导体 表面 电荷 采用 方程 如 下 


(7) 


4e0U0E,(U -Uo) 


0 


rHE.U E 一 四 


式 中 ，B, 为 导线 正 下 方 地 面 附近 的 标 称 场 强 ;7 为 
导线 半径 ， 歹 为 导线 高 度 ，U 为 导线 的 起 晕 电 压 。 
本 文 电荷 密度 迭代 公式 为 
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让 === 

求解 泊 松 方程 ， 得 

建立 模型 ， 给 出 到 各 单元 和 市 点 的 < 一 一 一 
单元 参数 电场 强度 


| i 


加 载 边 界 条 件 和 修改 怕 点 电 
初始 电荷 密度 茶 密 度 


| 


计算 各 单元 的 系 
数 和 矩阵 和 总 体 
合成 矩阵 


TE 


图 1 程序 计算 流程 


Fig.l The calculation process of the program 


所 Evra -Eb. 

Pi =p! | (9) 
式 中 , p? 为 布点 i 在 第 n 次 迭代 时 的 电荷 密度 ; 人 
为 修正 因子 ，E, 为 导线 表面 最 大 场 强 值 。 

本 文 值 选择 为 (Bw -5)/Ewse， 这 样 使 程序 在 
运行 时 根据 结果 能 够 自行 调整 表面 电荷 变速 度 。 在 
表面 误差 较 大 时 能 快速 向 收敛 发 展 ， 同 时 在 接近 收 
敛 值 时 ， 使 修正 变化 减 慢 ， 尽 可 能 提高 精度 。 

在 每 次 迭代 进行 修正 之 后 ， 再 次 进行 下 一 次 迭 
代 计 算 ， 在 满足 收敛 条 件 得 到 数值 解 后 ， 推 出 授 代 。 


4 ”仿真 结果 分 析 


4.1 离子 流 场 作用 下 空间 电荷 分 布 计算 

本 文 利 用 Matlab 软件 自 编程 求解 。 采 用 有 限 元 
法 计算 输电 线 附近 空间 电荷 之 后 ， 得 到 输电 线 附 近 
空间 电荷 的 分 布 。 

模型 的 原始 参数 为 ， 导线 半径 为 0.02m， 距 地 
高 度 为 20m; 起 始 电 晕 场 强 为 35kV/m， 起 始 电 晕 
电压 329KVG。 

本 文 分 别 计算 了 加 载 沿 x 轴 正 方向 2m/s、4m/s、 
6m/s 和 8m/s 的 风速 ， 计 算 不 同 风速 影响 下 输电 线 
附近 空间 电荷 密度 ， 结 果 如 图 2 ~ 图 5 所 示 。 从 图 
中 结果 可 以 发 现 水 平 风 速 从 0m/s 逐渐 增 大 到 4m/s 
时 ， 风 速 对 导线 附近 空间 电荷 分 布 存在 较 明 显 的 影 
响 ， 电 荷 密度 的 最 大 值 与 风速 呈正 比例 关系 ， 随 着 
风速 的 增加 ， 整 个 空间 电荷 密度 的 最 大 值 也 随 之 增 
加 ， 整 体 分 布 也 明显 向 风 的 方向 偏 移 。 

4.2 空间 电荷 对 绝缘 子 沿 面 电 场 参数 影响 
超 高 压 、 特 高 压 直 流 输电 系统 中 ， 由 于 电压 等 
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图 2 风速 2m/s 时 空间 电荷 密度 分 布 


Fig.2 The density distribution of the space charge in 2m/s’s 
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图 3 风速 4m/s 时 空间 电荷 密度 分 布 


Fig.3 The density distribution of the space charge in 4m/s’s 


wind speed 


级 高 ， 很 多 环境 因素 都 有 可 能 对 绝缘 子 的 电气 性 能 
造成 破坏 ， 导 致 绝缘 子 无 法 正常 工作 ， 因 此 对 绝缘 
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4 风速 6m/s 时 空间 电荷 密度 分 布 
Fig.4 The density distribution ofthe space charge in 6m/s’s 
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5 风速 3m/s 时 空间 电荷 密度 分 布 
Fig.5 The density distribution of the space charge in 8nmy/s’s 
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的 ， 本 章 计算 了 在 空间 电荷 存在 的 情况 下 绝缘 子 的 
电场 与 电位 分 布 ， 并 与 交流 线路 运行 状态 下 的 绝缘 


子 表面 电场 及 电位 进行 深入 的 研究 与 分 析 是 有 必要 。” 子 进行 了 对 比 。 模 型 具体 参数 见 下 表 。 
表 “绝缘子 主要 尺寸 和 参数 
Tab. The main dimension and parameter of the insulator 
盘 径 /mm 伞 盘 数 
额定 机 械 负荷 /kN A 蕊 棒 直 径 /mm 泄漏 距离 /mm 绝缘 长 度 /mm 单 边 金属 端 长 /mm 
160 144 ul0 24 54 54 12 786 4240 190 
; 电场 强度 /(V/m) 
在 无 空间 电荷 情况 下 电场 与 电位 分 布 由 图 6 和 i 时 ggrs 
加 BD 41/388 
图 7 所 示 。 加 945133 
> 0.104E+07 
在 存在 空间 电荷 条 件 下 ， 绝 缘 子 的 电场 及 电位 虽 拉夫 
分 布 云图 如 图 8 和 图 9 所 示 。 空 间 电位 由 输电 线 的 三 E+ 
电位 及 空间 电荷 共同 作用 决定 ， 可 以 看 出 由 于 空间 是 232 
电荷 的 作用 ， 整 个 空间 的 电位 有 明显 提高 晶 周二 1 
在 后 处 理 中 ， 将 绝缘 子 各 伞 群 边缘 点 的 电场 值 晶 和沙 


及 电位 值 提取 出 来 ， 图 10 为 包含 空间 电荷 及 不 包 
含 空间 电荷 的 电场 对 比 曲 线 ， 图 11 为 电位 对 比 曲 
线 图 。 


6 无 空间 电荷 电场 分 布 


Fig.6 The electric field distribution without space charge 
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7 无 空间 电荷 电位 分 布 
Fig.7 The potential distribution without space charge 
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图 8 存在 空间 电荷 电场 分 布 
Fig.8 The electric field distribution with space charge 
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图 9 存在 空间 电荷 电位 分 布 
Fig.9 The potential distribution with space charge 


由 图 10 可 知 ， 空 间 电荷 的 存在 使 伞 群 的 电位 值 
明显 升 高 。 同 时 ， 对 比 两 种 状态 下 的 金 裙 电 场 强度 
分 布 ， 可 以 看 出 考虑 空间 电荷 情况 下 ， 高 压 端的 电 
场 强度 明显 提高 ， 而 绝缘 子 低压 端的 电场 强度 有 一 
定 的 下 降 。 空 间 电 痪 的 存在 对 绝缘 子 表面 的 电位 和 
电场 强度 有 很 大 的 影响 。 产 生 这 种 现象 的 原因 是 导 
体 的 电 晤 层 具 相同 极 性 的 电荷 ， 在 空间 的 电位 有 空 
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10 ”绝缘 子 沿 面 电 位 分 布 对 比 


Fig.10 The contrase of the potential distribution along the 


surface of the insulator 


，。 不 存在 空间 电荷 
10| ， -存在 空间 电荷 


电场 强度 /( x 10 V/m) 


Xm 
图 11 绝缘 子 沿 面 电场 分 布 对 比 


Fig.11 The contrase of the electric field distribution along 


the surface of the insulator 


间 电 荷 和 ， 所 以 在 空间 上 对 绝缘 子 的 电动 势 有 一 定 
加 强 。 


本 文 仿真 分 析 了 500kV 特 高 压 直 流 输电 时 导线 
周围 空间 电荷 分 布 对 绝缘 子 沿 面 电场 分 布 的 影响 ， 
并 得 到 以 下 结论 : 

(1) 针对 土 500kV 单 极 线路 ,分析 了 风速 对 导 
线 附 近 空 间 电荷 分 布 的 影响 ， 风 速 会 使 空间 电荷 密 
度 最 大 值 升 高 并 向 风向 偏 移 ， 给 出 了 空间 电荷 分 布 
云图 。 

(2) 研究 了 存在 空间 电荷 情况 下 线路 上 熙 式 绝 
缘 子 的 沿 面 电位 及 场 强 分 布 ， 并 进行 了 分 析 。 绝 缘 
子 高 压 端的 电场 强度 小 于 静电 场 的 电场 强度 ， 绝 缘 
子 低 压 端的 电场 强度 大 于 静电 场 的 电场 强度 ， 大 约 
增加 了 8% 。 
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